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Abstract 

This paper is devoted to study the accelerated expansion of the 
universe by exploring the Brans-Dicke parameter in different eras. For 
this purpose, we take FRW universe model with viscous fluid (without 
potential) and Bianchi type I universe model with barotropic fluid 
(with and without potential). We evaluate deceleration parameter 
as well as Brans-Dicke parameter to explore cosmic acceleration. It 
is concluded that accelerated expansion of the universe can also be 
achieved for higher values of the Brans-Dicke parameter in some cases. 
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1 Introduction 



The accelerated expansion of the observable universe is one of the most con- 
spicuous and recent achievement in modern cosmology. This expansion with 
positive cosmic acceleration has been confirmed by many astronomical ex- 
periments such as Supernova (la) l ' 2 \ WMAP 3 \ SDSS 4 \ galactic cluster 
emission of X-rays 5 \ large scale structure 6 -*, weak lensing 7 > etc. These 
results lead to the conclusion that our universe is spatially flat. 
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The positive cosmic acceleration of the universe has been motivated by 
a mysterious exotic matter having large negative pressure, known as dark 
energy. Although, General Relativity (GR) is an excellent theory to explain 
the gravitational effects but it is unable to describe the present cosmic ac- 
celeration and the reality of dark energy. In order to explain the nature of 
this mysterious finding, various models including Chaplygin gas, phantom, 
quintessence, cosmological constant etc. have been constructed 8 ' 9 \ However, 
none of these models is very successful. 

The exploration of scalar-tensor theories of gravity as modified theories of 
gravity has received much attention due to their vast implications in cosmol- 
ogy 10 ~ 14 ). Brans-Dicke (BD) theory of gravity, a special case of scalar-tensor 
theories, is one of the most viable theories for this purpose. It is the general 
deformation of GR satisfying weak equivalence principle, in which gravity ef- 
fects are mediated by the metric tensor and the scalar field 15 \ This provides 
a direct coupling of the scalar field to geometry. Brans-Dicke theory is com- 
patible with both Mach's principle 16 ) and Dirac's large number hypothesis 
17 \ One of the salient features of this theory is that the gravitational cou- 
pling constant, being the inverse of spacetime scalar field, varies with time. 
In order to fulfill the solar system experiment constraints, the value of the 
generic dimensionless BD parameter oj should be very large, i.e., u > 40, 000 

18,19) 

Brans-Dicke theory is a successful theory which can tackle many out- 
standing cosmological problems like inflation, quintessence, late time behav- 
ior of the universe, coincidence problem, cosmic acceleration u ) etc. There 
are different versions of BD theory available in literature 20 ' 21 ). Singh and 
Rai 22 ) investigated various BD cosmological models and showed that the 
Bianchi models are very effective in explaining the evolution of the universe 
for perfect fluid. Bermann 23 ) discussed different models of the universe with 
constant deceleration parameter based on the variation law of Hubble param- 
eter. Bertolami and Martins n ) found that the accelerated expansion of the 
universe could be obtained with large \u>\ and potential 2 without consider- 
ing the positive energy condition. Sen and Sen 24 ) showed that the dissipative 
pressure could support the late time accelerated expansion of the universe. 
Banerjee and Pavon 12 ^ found that the present accelerated expansion could 
be obtained without restoring a cosmological constant or quintessence matter 
for FRW model. 

Sahoo and Singh 25 ) explored the observed accelerated expansion of the 
present universe in this theory for FRW model. The same authors 13 - ) found 
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exact solutions in different eras of the universe and discussed the possibil- 
ities for obtaining cosmic acceleration, inflation and deceleration for these 
solutions. Sen and Seshadri 26 ) investigated the role of positive power law 
potential on the accelerated expansion of the universe. They concluded that 
self-interacting potential can derive the accelerated expansion in the perfect 
fluid background with small negative values of BD parameter. Reddy and 
Rao 27 - ) found axially symmetric perfect fluid cosmological model in this the- 
ory. In order to investigate the present accelerated expansion of the universe 
and different stages of the cosmic evolution, a lot of work has been done using 
Bianchi models in GR and scalar tensor theories 28 ~ 32 ). In a recent paper, 
Chakraborty and Debnath 14 ) have investigated cosmic acceleration in this 
theory for FRW model. They have shown that the accelerated expansion of 
the universe with higher values of u can be achieved only for closed model. 

In this paper, we explore the role of BD parameter on the cosmic accel- 
eration by using spatially flat models in the presence of different fluids. The 
paper is organized as follows. In the next section, we formulate the field 
equations of generalized BD theory with self-interacting potential. Section 3 
provides the field equations for FRW model in the presence of a viscous fluid. 
Here we discuss models for both constant as well as varying bulk viscosity 
coefficient. In section 4, we formulate the field equations in the presence of 
the barotropic fluid for the Bianchi type I (BI) universe model. This section 
explores all possible choices of the BD parameter u and the self-interacting 
potential V(<f)). Section 5 investigate the observational limit of gravitational 
constant for the constructed models. Finally, we discuss the results in the 
last section. 



2 Brans-Dicke Field Equations 

Brans and Dicke 15 - ) proposed a scalar-tensor theory known as Brans-Dicke 
theory of gravity that was based on the pioneering work of Jordan. A mod- 
ified version of this theory is the generalized BD theory in which the BD 
parameter no longer remains a constant rather, it turns out to be a function 
of the scalar field. The action for generalized BD theory with self-interacting 
potential in Jordan frame 20 ' 21 ) is given by 

S = [ d 4 xV=g[<pR - ^P-<P' a (P, a ~ V{4>) +L m ]; a = 0, 1, 2, 3, (1) 



3 



where BD parameter u(4>) is the modified form of the original BD parameter 
oo, V(4>) denotes the self-interacting potential and L m represents the matter 
part of the Lagrangian. Here we have taken 87rGo — c — 1. Taking variation 
of this action with respect to the metric tensor g^ v and the scalar field, we 
obtain the following BD field equations 14 ** 



2<p 



9r + ^> (2) 



□0 = — L.w)-^] 

r 3 + 2w(0) 3 + 2w(0) [ W) V d(j> - 



dui((j>) 
d</> 



Au^, (3) 

where T = g^T^ denotes trace of the energy-momentum tensor and □ 
represents the d'Alembertian operator. Equation ([3]) is called wave equation 
for the scalar field. Notice that BD theory is reducible to GR if oo — > oo and 
the scalar field becomes a constant 33 ^. However, it is not true in general. 
In papers 30 > 34 ) ; it has been pointed out that BD theory does not always go 
over to GR in the limit oo — > oo for the case of exact solutions. In this limit, 
GR could be recovered only if the trace of the energy-momentum tensor 
describing all fields other than BD scalar field does not vanish, i.e., ^ 
34-37) _ p or rp(j ri ) _ gj) solutions do not correspond to respective GR 

solutions. The Palatini metric f(R) gravity and the metric f(R) gravity are 
obtained by substituting oo = —3/2 and oo = respectively 38 ^. 



3 Cosmic Acceleration and FRW Model 

In this section, we investigate cosmic acceleration by exploring the BD pa- 
rameter. For this purpose, we consider FRW model with viscous fluid. In 
particular, we discuss two cases according to bulk viscosity is constant or 
variable. The line element for FRW model is given by 

dv^ 

ds 2 = dt 2 - a 2 (t)[ - + r 2 (d6 2 + sin 2 0d<t>% (4) 

1 rCV 

where a(t) is the scale factor and k = —1, 0, +1 indicate open, flat and closed 
universe model respectively. We assume that the universe is filled with a 
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viscous fluid given by 



Tfj, v = (p + P e ff)u^U v - Peffg^u, (5) 

where p is the energy density, is the four-vector velocity satisfying the 
relation u^u" = 1 and P e ff represents the effective pressure defined by 

Peff = Pi + P v is- 

Here Pi denotes the isotropic pressure and P vis represents the pressure due 
to viscosity. The bulk viscous pressure is defined by Eckart's expression 
in terms of fluid expansion scalar and is given by P V i S = — 39 \ where 
£ = £(*)P) represents the bulk viscosity coefficient. For FRW model, the 
viscous pressure is found to be P V i S = — — and hence the effective pressure 
becomes 

P eff = Pi~ 3£#, (6) 

where H = - denotes the Hubble parameter. The corresponding field equa- 
tions turn out to be 

a 2 + k a<p u <fi 2 p 

~^ 2 ~ + ~^4>~ 6 02 = 30' ( > 



2a a 2 + k u 2 2d0 , -P t 
a a 2 2 2 a0 



where dot denotes the derivative with respect to time. The corresponding 
wave equation becomes 

+ 3tf0=^^-^. (9) 

Y r 2w + 3 2w + 3 v ' 

Here we have taken V(<f)) = 0. 

Equation of state provides a relation between isotropic pressure and en- 
ergy density and is given by 

P = IP, (10) 

where 7 is the equation of state (EoS) parameter. The values of 7 = 
— 1, 0, 1/3, 1 represent vacuum dominated, dust, radiation dominated era 
and massless scalar field respectively. The continuity equation for the viscous 
fluid (JSJ) can be written as 

p + 3H(p + P eff ) = 0. (11) 
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One can assume the standard expression for bulk viscosity, i.e., £ = £oP n > 
where n is a non- negative constant and £o > 0. Different possible values 
of n are available in literature 40_43 ) ; out of which two choices n = 1 and 
n = 3/2 correspond to the radiative and string dominated fluids respectively. 
However, more realistic models can be obtained for < n < 1/2. Here 
we would like to evaluate p by solving the continuity equation ffTTj) for the 
following two cases: 

• Constant bulk viscosity, i.e., £ = £o (f° r n — 0). 

• Variable bulk viscosity, i.e., £ = £(t,p) with n = 1/2, 1. 
In both cases, we choose k = 0, i.e, flat FRW model. 



3.1 Constant Bulk Viscosity Coefficient 

The energy conservation equation ( jTTj) . in terms of constant bulk viscosity, 
can be written as 

p + 3H(l + 1 )p = 9^H 2 , (12) 

where we have used the EoS given by Eq.f JTOj) . We assume that the scale 
factor a(t) has the form of expanding solution (power law form) 

a(t) = a t a , a > 0, (13) 

where is the present value of the scale factor. The deceleration parameter 
is given by 

9= + (14) 

Notice that the deceleration parameter q suggests a > 1 for cosmic acceler- 
ation. Equation (1121) leads to 

pit) = ^! r i + Po a 3 ^h~ 3 ^K (15) 

V ' [-1 + 3a(l + 7)] H V ' 

The scalar field can be found from Eq.Q by taking cj(t) = u (constant) as 
follows 

(1 - 3 7 )p a 3(1+7) ^ 2 " 3a(1+7)1 , 3£ (l-4a)t 



(3 + 2w )(l - 3«7)[2 - 3a(l + 7)] (3 + 2u )[l - 3a(l + 7)] 
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This equation suggests that the scalar field can be taken in the power law 
form when the scale factor is given in the expanding form. 

We would like to discuss the time dependent BD parameter u which 
satisfies the field equations as well as wave equation. For this purpose, we 
assume a simple form of power law for the scalar field 

<f)(t) = fat?, (16) 

where 0o is the present value of the scalar field and 3 is any non-zero constant. 
The field equation (jTJ) can be re-arranged in the following form 



a j Mt)P | P (u) 
a 20 V 1202 +30' V) 

where fi(t) = 2u(<f>(t)) + 3. Using Eqs.(fl3]) and fifty in Eq. (117)1. it follows 
that 

= 3^[^ - (18) 
The comparison of Eqs. f|T5|) and f fT8|) yields 

Q( t ) = —(2a + B) 2 t {l ~ p) - 1 __ 

iJ ^[3a(l + 7 )-l] o/32 

(19) 

The corresponding expression for will become 

" (0(t)) "" (t) = W(l + 7)-l] «P ' (20) 

Here we have considered the time dependent terms only. 

In order to check the consistency of these solutions with the wave equa- 
tion, we substitute these values in (J5J). This leads to the following two con- 
sistency relations 

B[A(8 - 1 + 3a) + B(l - 3 7 ) + 2 [2 - 3a(l + 7) - P]] = 0, (21) 
4a(/3 - 1 + 3a) + a(l - ^)B + [3a (1 + 7) - 1]B + 2a(l - B) = 0. 

(22) 

Equation (j2Tp implies that either /3 = or /3 = —2a while Eq. fT22l is satisfied 
for either B = —2a or a = 1/6. For cosmic acceleration a = 1/6 is not an 
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interesting value and so we ignore it. When = 0, the BD parameter yields 
uj(t) — > — oo and the scalar field becomes a constant, i.e., = O . This 
leads to GR, so it is not the interesting case. For = —2a, u takes the form 

nt n +(l+2«) -3(1+7), [2-0,(1+37)] 

U 20 o [3a(l + 7)-l] 20 o a 2 ' 1 } 

The power law expression for the scalar field turns out to be 

<f>(t) = U~ 2a . 

In the following, we evaluate the BD parameter at different epochs of the 
universe. 

For vacuum dominated era (7 = —1), the BD parameter is 

u{t) = ^t^--^- 2 t^\ (24) 
200 200"^ 

while in the radiation dominated era (7 = 1/3), it becomes 

t " W "20 o (l-4«) t 20 o «2 ■ {2b) 

The BD parameter in the matter dominated era or the dust case (7 = 0) 
takes the form 



In the massless scalar field era (7 = 1), this turns out to be 

^ j ~20 o (l-6a) 20 o a 2 • W 

Finally, the BD parameter for the present time, t = to, can be calculated 
from dust case, i.e., matter with negligible pressure. Equation ( 126]) leads to 
the present value of the BD parameter w given by 

U ° ~ "2(3q-1) " 2rf- (28) 
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Figure 1: Plots show the graph of u(t) versus cosmic time t for (a) 7 = — 1 
and (b) 7 = 1/3 with a = 1.1 and varying values of £0 as follows: blue, 
£0 = 0.0001; red, £ = 0.15; pink, £ = 0.2; green, £ = 0.38. 



Here we have normalized the constants, i.e., (f) = a = t = 1, p = and 
a > 1 which is consistent with Eq.f ri7|) . The minimum value of ooq is 



- _^ _ I 
U ° 4 2' 

Clearly the minimum value of ujq depends on the value of constant bulk 
viscosity coefficient £o- In the BD theory, the gravitational coupling constant 
and the scalar field density should be positive in the present universe which 
can be achieved for u; > —3/2 12 \ In our case, the bulk viscosity coefficient 
must have £0 < 4/9 with a > 1 for the consistency purpose. The present 
observational range for the deceleration parameter is —1 < go < 1>2 ^ which 
restricts a > 1. The more general form of the model for the present universe 
can be obtained by taking a — 1 + e, e>0 (for small values of e) given by 

u f t ) = -4^ - -J = o r 2(1+e) , a(t) = a t (1+£) . 

W 2(2 + 3e) 2(1 + e) 2 ' r ro ' w 



Now we discuss the BD parameter for vacuum as well as matter domi- 
nated eras. In the vacuum dominated era, the graphs indicate that u(t) is a 
decreasing function starting from zero for < £0 < 0.11. For £0 > 0.11, the 
graphs correspond to increasing function but after some particular points, 
they again become decreasing function as shown in Figure 1. Thus, for this 
range of constant viscosity £ with a > 1, it is possible to achieve the cosmic 
acceleration with positive values of u(t). In all other eras of the universe, u(t) 
is a decreasing function of time with smaller negative values. For < a < 1, 
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in the radiation dominated era, u(t) is a decreasing function and the universe 
undergoes to decelerated expansion. Thus, £ plays the role to control the 
time dependence of u(t). In the radiation, matter dominated eras and mass- 
less scalar field, the BD parameter approaches to — oo for a = 1/6, 1/3 and 
1/4. For the cosmic acceleration, we must have a > 1, hence these values 
are not interesting. 



3.2 Variable Bulk Viscosity Coefficient 

For the sake of simplicity, we take n = 1/2, i.e., £(t, p) = iop 1 ^ 2 ^). Using 
this value of bulk viscosity coefficient along with Eq. (I13p in ( II ip , it follows 
that 

r 3a N 9£ ct 2 



p(t) + [-j- (l+ 7 ) f . 



-} P (t) 1/2 = o. 



This yields the following solution 



P(t) 



t[3a(l + 7 ) -2] 



+ Pot 



-3a(l+7)/2 



(29) 



where po is an integration constant. Comparing this equation with Eq.(IT8l). 
we obtain 



2a + /3 



At 2 



t[3a(l + 7) -2] 



(3 J ^o/? 2 ^ 
The corresponding BD parameter will become 

2t 2 T 9£ a 2 



+ Pot 



-3a(l+ 7 )/2 



<j) P 2 t? [t[3a(l + 7) - 2] 



(30) 



For consistency of this solution with the wave equation, we substitute all 
these values in the wave equation (jUJ) which leads to 

a = 1/3, (3 = —2a, 7 = 1. 

For (3 = —2a, we obtain 



u{t) 



t 2(l+a) 

' 20 o « 2 



9eo« 2 



t(3a(l + 7) -2) 



+ Pot 



-3a(l+7)/2 



(31) 
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Figure 2: Plots show the graphs for u(t) versus cosmic time t for a = 1.1, (a) 
7 = and (b) 7 = 1/3 with varying values of £o as follows: red, £o = 0.0001; 
green, £ = 0.38; blue, £o = 0.01; pink, £ = 0.1. 



The choice a = 1/3 is not feasible for obtaining cosmic acceleration while 
7 = 1 corresponds to massless scalar field which is discussed below. Now we 
evaluate BD parameter for different eras. 

In the vacuum dominated era, the BD parameter is 



u{t) 



2<p a 2 



2t 



Po 



while for the radiation dominated era, it turns out to be 



t 2(l+a) 

' 2<p a 2 



9£o« 2 



t(4a - 2) 



+ Pot 



-2a 



The BD parameter in the matter dominated era is 



t 2(l+a) 



-i 2 



t(3a - 2] 



+ Pot 



-3a/2 



2(j) a 2 

For the massless scalar field era, this is given by 

£2(l+a) r r.^ 9 n2 



U{t) 



20 o « 2 



t(6at - 2) 



+ Pot 



-3a 



(32) 



(33) 



(34) 



(35) 



The expressions for u(t) correspond to decreasing function as — t 2 for in- 
creasing values of viscosity coefficient £o and — 1 < 7 < 1 except for vacuum 
dominated era. This gives rise to accelerated expansion of the universe for 
a > 1 as shown in Figures 2 and 3. For a = 1/2, 3/2 and a = 1/3, the BD 
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Figure 3: Plots shows the graphs for u(t) versus cosmic time t for a = 1.1 
and 7 = — 1 with varying values of £o as follows: red, £o = 0.0001; green, 
£o = 0.38; blue, £ = 0.01; pink, £ = 0.1. 

parameter approaches to — oo. Here in the matter dominated era, a = 3/2 
lies in the range a > 1 allowed for the accelerated expansion of the universe. 

Now we discuss the radiative fluid case (n = 1). Here we take £(t, p) = 
£op(t). Consequently, the continuity equation yields 

p(t)=Pot- 3Q(1+ ^exp(^|^). (36) 

The BD parameter u(t) turns out to be 

U f t) = 2p exp(^) ^ 2 _ 3a(1+Tl _ g1 
4>oP 2 

Here (3 = and /3 = —2a are the corresponding consistency relations. The 
choice of /3 = provides no interesting insights while = —2a leads to the 
following expression 

u{t ) = _^P(f5zl t [2-a(l + 3 7 )]_ (37) 

20 o « 

For the radiation dominated era, the BD parameter takes the form 

= _ poexp(y) f2 „_ e ,, 
2(p a z 

We see that the coefficient of viscosity appears only in the exponential func- 
tion. Thus, in the radiation dominated era, for small values of £o an d 
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Figure 4: Plots show the graphs for u(t) versus cosmic time t for (a) a = 1.1 
and (b) a = 0.5 with varying values of £o as follows: green, £ = 0.0001; red, 
£o = 1-2; blue, £ = 0.09; pink, £ = 0.38. 

a > 1, exp(— 9£ «A) — > 1, providing small negative values of u(t) as shown 
in Figure 4. If £ — > oo with a > 1, then exp(— 9£oa/t) — > which means 
cj(t) — > 0. Thus this model may correspond to that of the metric f{R) grav- 
ity. However, it is not physically possible. Also, in this case, for < a < 1, 
the values of u(t) are constrained within the range —3/2 < u(<j)) < which 
shows that the universe undergoes to the decelerated phase. 

4 Cosmic Acceleration with Barotropic fluid 
and Bianchi I Universe Model 

Here we investigate expansion of the universe by using LRS Bianchi type / 
model in the barotropic fluid background. The line element of Bianchi type 
I universe model is described by 44 -* 

ds 2 = dt 2 - A 2 (t)dx 2 - B 2 (t)(dy 2 + dz 2 ), (39) 

where A and B are the scale factors. This model has one transverse direction 
x and two equivalent longitudinal directions y and z. Assume that matter 
contents of the universe are described by the perfect fluid given by 

T,,,, = + PWu u - P 9uv . (40) 
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The corresponding field equations (j2J) and (EJ) can be written as 

AB + B^~ ( j ) + 2 02+ 20 M + V0' 1 J 

2^ + 5! = _ ^M^! _ 2 ^ - * + f42l 

5 5 2 2 2 5 20 ' 1 ' 

^,^ + M = _^_ £ _ i + YM. _ r d + i\i (4*\ 

B A AB 2 2 20 [ A B } <f)' 1 1 



The wave equation is 



$ +( A +2 *)i = p~ 3P J 2V ^-* « i ^_ (44) 

^ V A 2w(0) + 3 2^(0) + 3 2w(0) + 3' V ; 

For this model, the average scale factor and the mean Hubble parameter are 

a\t) = AB\ H(t)= 1 -(j + 2^). 

The energy conservation equation for energy-momentum tensor given in 
Eq.gQD will be 

P + 4 + 2§)(p + P) = 0. (45) 

We assume that the universe is filled with barotropic fluid. The barotropic 
EoS 14 - ) is given by 

P = ip; (-1<7<1)- 

The expansion scalar for Bianchi type I model is given by 
while the shear scalar is 

1 A B 
a= ^ { A-B ] - 

It is given 45 - ) that for spatially homogeneous metric, the normal congruence 
to homogeneous expansion yields the ratio | as constant i.e., "expansion 
scalar 9 is proportional to shear scalar a". This physical condition leads to 
the following relation between the scale factors 

A = B m , (46) 
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where m ^ 1 is any positive constant (for m = 1, it reduces to flat FRW 
model). In literature 44_49 ) ; this condition has been widely used to find exact 
cosmological models. Using this assumption in Eq. (I45[) . it follows that 

p+(l+ 7 )(m + 2)|p(t) = 

which yields 

p(t) = p £- (1+7)(m+2) . (47) 
Now we discuss the various possible choices for oj(4>) and V(0). 

4.1 Model without potential V{<j>) = 

We take the following two cases according to u> is constant and u = u>{4>). 

4.1.1 Case (i) 

First we take BD parameter as a constant, i.e., u)(<f>) = u . For the solution 
of the field equations, we consider the power law as follows 

B(t) = b t a , a>0. (48) 

Using Eqs. (l46p . (I48p and the mean Hubble parameter H, the deceleration 
parameter can be written as 

q=-[l- 



a(m + 2) J ' 

Notice that q < 0, q = and q > indicate an accelerated expansion, 
uniform expansion and the decelerating phase of the universe respectively. 
Thus, for accelerated expansion of the universe, we must have the following 
condition on a 

3 

a > ; m^l. (49) 

[m + 2) 

Substituting Eqs. (l4"S]) and (14"8"]1 in (B3J, the scalar filed becomes 

(1 - 3 7 )p &-( m+2 )( 1+ ^tP-"(l+7)( m+ 2)] 

0U (3 + 2wo)[l-a7(^ + 2)][2-a(l + 7)(m + 2)]" 1 j 
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Figure 5: Plot of loq versus a for (a) m = 3/2, (b) m = 4/5 and 7 = — 1. 
The corresponding ranges for a are a > 6/7 and a > 15/14. 



(a) For y=-4/5 and m=3/2 



u (b) For y=-4/5 and m=4/5 

■1 ■ 

2 



Figure 6: Plot of coo versus a for (a) m = 3/2, (b) m = 4/5 and 7 = —4/5. 
The corresponding ranges for a are a > 6/7 and a > 15/14. 



The BD parameter is obtained from the field equations ( )4T|) - (l43|) as 

1 (m + 3)a(q - 1) 2 [(m 2 + 1) + 2 7 (2m + 1)] 

"(T^) [ [2-a(m + 2)(l + 7 )] 2 + " [2 - a(m + 2)(1 + 7 )] 2 
q[m + 3 + 27(m + 2)] [1 - q(m + 2)(1 + 7)] 
+ [2-q(m + 2)(l + 7 )] + [2-q(m + 2)(l + 7 )] J - ( ' 



For massless scalar field 7 = 1, we have u — > —00, which leads to GR. 
We have seen that the BD parameter depends upon the parameters a, 7 
and m. These parameters are constrained using some physical conditions. 
The possible ranges for m are < m < 1 and m > 1 and 7 is allowed 
for — 1 < 7 < 1. The deceleration parameter constraints a such that a > 
3/(m + 2). By taking different possible choices for these parameters, it can 
be seen that BD parameter takes small negative values as well as positive 
values for —1 < 7 < as shown in Figures 5-7. This gives rise to cosmic 



16 



w (a) For y=-9/10 and m=3/2 " (b) For y=-9/10 and m=4/5 




a 



Figure 7: Plot of ujq versus a for (a) m = 3/2, (b) m = 4/5 and 7 = —9/10. 
The corresponding ranges for a are a > 6/7 and a > 15/14. 



acceleration for this range of 7. We would like to mention here that for ceratin 
ranges of a allowed for cosmic acceleration and — 1 < 7 < 0, uq can take 
larger values which would be compatible with the solar system experiment 
constraints. 

Solving Eq. (|5ip for a, we obtain the following quadratic equation 

a 2 [(m + 2) 2 (1 + 7 ) 2 [( 7 - l)w - 2] - (m + 3) - [m 2 + 1 + 2 7 (2m + 1)] 
+ (m + 2)(1 + 7 )[m + 3 + 2 7 (m + 2)]] + a[{m + 2)(1 + 7)[~4(7 - l)w 
+6](m + 3) - 2[m + 3 + 2 7 (m + 2)]] + 4[w (7 - 1) - 1] = (52) 

which provides two roots. These values for m — 1/2 and 7 = (present 
universe) are given by 

23/2 + 10u;o± ^-15/4-6^0 

a = —, . (53) 

17 + 25/2u v ; 

Since — 2 < u Q < — 3/2 is the observed range for cosmic acceleration, so the 
choice of Uq = —5/3 leads to following values of a 

ax = 16/23, a 2 = 2. 

Here a± gives q > 0, hence we leave it while a 2 = 2 yields g < leading 
to accelerating expansion. Also, it yields <p(t) = t~ 3 which provides positive 
coupling constant. Since in our case, 4>(t) is decreasing more rapidly as 
compared to <f){t) = t~ 2 n > and <f>(t) = t~ 5 ^ 2 l2 \ therefore it corresponds to 
greater rate of accelerated expansion of the universe. 
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4.1.2 Case (ii) 

In this case, the BD parameter is not constant rather it is a function of 0. 
Using Eqs.((l3D, flU, (|Hj) - (|13)» and flU, the BD parameter can be written 

as 

1 (3m -m 2 -2) 2 (m + 3)a (m + l)q/3 2 
W (<W = ^1 2 2 2 ^ ^ 

1 ^ lOWl , \ -c(m+2)(l+ 7 )-/3+2 -"("»+2)(l+ 7 )-2 

^[po& (m+2)(1+7) (1+7)0 " 0o " ]• (54) 

Substituting this value in Eq.f l44l) . we obtain the following consistency rela- 
tion 

fl = _(m + 2)a(l + T ), m#1 (M) 

This shows that /3 remains negative for all < m < 1, m > 1, a > 3/(m+2) 
and — 1 < 7 < 1. The consistency of this solution with the dynamical 
equations, i.e., each term in the dynamical equations should have the same 
time dependence, results in another constraint given by 

/3 = 2-a(m + 2)(l + 7 ). 



Using this value of (3 in Eq. fl55|) . it can be seen that the parameter f3 is 
restricted to —2. Now we discuss the BD parameter and cosmic acceleration 
in different phases of the universe by using this value of 0. The expressions for 
BD parameter in matter and radiation dominated eras with (3 = — 2, a > 6/5 
and m = 1/2 turn out to be 

-(0) = ^hf ! + ^-6]-^- 2+ n 

M) = + 6 - -6 { - 2+ ^\ 

^ 4 l 8 4 J 3^ 

By taking different choices for these parameters, we see that for all phases 
of the universe, the BD parameter u(<f)) has small negative values and lies 
within the range u < —3/2 as shown in Figures 8 and 9 which corresponds 
to accelerated expansion of universe. This result is in agreement with 14 -* for 
spatially flat model. 

4.2 Model with potential V(<j>) ^ 

Again, we discuss two cases depending upon the value of BD parameter u. 
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Figure 8: Plot of oj(<p) versus for m — 1/2 and (a) 7 = 1/3 and (b) 7 = 
with a > 6/5 as follows: red, a = 1.3; blue, a = 2; pink, a = 1.5. 
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Figure 9: Plot of oj{<p) versus <fi for m = 1/2 and (a) 7 = —1/3 and (b) 
7 = —1/2 with a > 6/5 as follows: red, a = 1.3; blue, a = 2; pink, a = 1.5. 
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4.2.1 Case (i) 

First we discuss the case of constant BD parameter, i.e., lo{4>) = <^o- Further, 
we consider the power law form of the scalar field in terms of scale factor 
Bit) 

4> = 4> B a ; a>0. (56) 
Using this value of 4> in the field equations (HT|) - fH5|) . it follows that 

R R 2 

_ + A(-) = _(7R-«(l+7)(m+2) 

B B 2 



where 



(a 2 — 3a — 2 + 2m 2 + am — 2m + u a 2 ) (l + j)p 

J\ = —~ ~~ : , O 



(3a + 2m) 0o(3« + 2m) 

The expression for B(t) can be written as 



B(t) = x B{t) 2 ~ { ^ )(m+2) x [(3a + 2m)[(l + 7 )(m + 2) 

+ (a -2)] - 2a 2 + 6a + 4 -4m 2 - 2ma + 4m - 2w a 2 ]" 1/2 (57) 
which yields 

B(t) = A't^ a+{l+ ^ m+ ^\ (58) 

where 

A> = [^[a + (l + 7)(^ + 2)] 2 (l + 7) x [(3a + 2m)[(l + 7 )(m + 2) + (a 
200 

- 2)] -2a 2 + 6a + 4- Am 2 - 2ma + 4m - 2u a 2 ]- 1 ] 1/la+{1+ ^ {m+2)] . 

(59) 

The value of the scale factor A(t) can be obtained by using value of B(t) in 
Eq.flU]). 

The corresponding expression for the scalar field is given by 

0(t) = ot 2a/[a+(l+7)(m+2)] ; ^ 

where 6 Q = <poA' a . Equation ( 15 8 p yields the following constraint on a 

3a < -(l + 37)(m + 2). (61) 
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The deceleration parameter q turns out to be 



_ 3[g + (l + 7 )(m + 2)] 

q 2(m + 2) 

It can be easily seen that for all positive constant m (m ^ 1), a > 1 and 
— 1 < 7 < 1, the deceleration parameter remains negative, i.e., g < 0. Thus 
the universe is in the state of accelerated expansion. From the wave equation 
flSJ), the potential can be written as 

B' 

V &) = (i +7 ) (m+ 2)/^ ( 62 ) 



where B' is given by 



— a>9, 



2/(op) 

o 



^' = 7 T7 Nr , -r[4(l + 7)(m + 2)(w a 

(l + 7)(m + 2)[a + (l + 7 )(m + 2)] 2L V nK A0 

- m - 2) - 8(m + 2)aw + 16m + 8m 2 + 24 - 4ma - 8a]. 
4.2.2 Case (ii) 

Let us take the BD parameter as a function of the scalar field 0, i.e, u)(<f>). 
Consider the power law forms for the scalar field and the scale factor, given 
by (l4"8]) and (fl3"|) . Using the field equations (|4T|) - (143|) . the scalar potential 
takes the form 

VU) = ^V 1 ^[ (m2 + 5 ' 2 " + 6)a2 - ( m + 3 >% /?-g + ( ( 3m + 7 > ),»/?] 
- 2 (63) 



The BD parameter turns out to be the same as given by f[54|) . Substituting 
these values in Eq.(l44p. we obtain the following consistency relations 



a 

/3 = 0, p = -2 /3 = --(m + 3), (3 = l- a{m + 2). (64) 

For the consistency of this solution with the dynamical equation, it follows 
that 

/3 = 2-a(m + 2)(l + 7 ). 
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Figure 10: Plot (a) shows the graph of V(4>) versus scalar field for m = 1/2, 
a = 1.3 and (3 = —7a/ '4 as follows: red, 7 = 1/3; green, 7 = — 1; pink, 7 = 
and Plot (b) shows the graph of V((f>) versus scalar field <fi for m = 3/2, 
a = 1.3 and (3 = — 9a/4 as follows: blue, 7 = 1/3; green, 7 = —1; red, 7 = 0; 
pink, 7 = 1. 



Now we discuss the behavior of self-interacting potential V(4>) for these 
values of in different eras of the universe. The choice /3 = is not feasible, 
so we neglect it. For (3 = —2, the self-interacting potential can be written as 

V&) = + f 1 + - )a 2 -(^)«+6-(3m+7) a ]-(l- 7 )^ + ^ 1 ^/ 2 i 

where m 7^ 1 is a positive constant and a > 3/(m+2). For (3 = — ^f, m = 3/2 
and a > 6/7, we get 

V{(j>) = 0< 1+ ^[-3a 2 ] - (1 - 7 )0 14 ( 1+ ^ 9 . 
For m — 2, a > 3/4 and (3 = — the potential turns out to be 

V(0) = -(l- 7 )0 8(1+7)/5 , 

where — 1 < 7 < 1. The expression for self- interacting potential for radiation 
dominated era with (3 = 2(1 — a > 6/5 and m = 1/2 is given by 

,95a 2 13a 4 N -W3 

1/(0) = ( h -)0(i-W3) . 

72 4 3 ;V 

The self interacting potential for matter dominated era with /3 = 2 — a > 
6/5 and m = 1/2 takes the form 

V(0) = (1 - — )0(2~W2). 
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Figure 11: Plot (a) shows the graph of V(<p) versus scalar field for m = 1/2, 
7 = 1/3 and (3 = 1 — 5a/2 as follows: blue, a = 1.5; green, a = 2; pink, 
a = 1.3 and Plot (b) shows the graph of V(4>) versus scalar field for m = 3/2, 
7 = and = 1 — 7a/2 as follows: blue, a = 2; green, a = 1.3; red, a = 1.5. 



v (<*) For 0=-2 




Figure 12: Plot shows the graph of V(4>) versus scalar field for m = 1/2 and 
(3 = — 2 as follows: blue, a = 1.5; red, a = 1.3 (7 = 1/3); green, a = 1.5; 
pink, a = 1.3(7 = — 1); sky blue, a = 1.5; yellow, a = 1.3(7 — !)• 
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V axis ( a ) V axis (j^ 




Figure 13: Plot (a) shows the graph of V(4>) versus scalar field for m = 1/2, 
7 = 1/3 and (3 = 2(1 — 5a/3) as follows: blue, a = 1.5; green, a = 2; 
pink, a = 1.3 and Plot (b) shows the graph of V(<j>) versus scalar field for 
m = 1/2, 7 = and (3 = 2 — ha/2 as follows: blue, a = 1.5; green, a = 2; 
pink, a = 1.3. 



For the first three consistency relations for f3 given by Eq. (l64p . we see that 
V(4>) is a decreasing function starting from zero with the increasing values 
for <fr except for the choice (3 = —2. In this case, only 7 = — 1 and 7=1 
with a > 1.3 provide positive potential energy as for these ranges, they are 
increasing functions of <fi as shown in Figures 10-12. Figure 13(a) shows that 
V(4>) attains negative values starting from zero but with larger values for a, it 
is an increasing function with positive values. Figure 13(b) shows that V^(^) 
attains positive increasing values for a > 6/5. Therefore, we conclude that 
these cases provide positive potential energy as they are increasing function 
of 4> for particular values of a. 

5 Variation for the Newton's Gravitational 
Constant in GBDT 

A well-known fact about BD theory of gravity is that it provides very small 
variations for the gravitational constant. However, GBDT suggests various 
possibilities for variation of G. In GBDT, the expression for G is found 20 ) 
to be 

G(t) = 4 + 2u ^ 
U 0(3 + 2^(0))' 
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The present rate of variation for gravitational constant is given by 





2(^o) f65) 
(3 + 2w )(4 + 2u;o)' 

Here subscript indicates the present values of the corresponding parameters. 
Using Eq. ([2"3]) . f3 = —2a, a > 1, £ = 0.0001 and the estimated age of the 
universe t = 14 ± 2 Gyrs, we obtain the rate of variation of (G/G)q to be 
1.5714 x 10~ 19 yrs. It lies clearly within the allowed range of variation of G 
for cosmic acceleration, that is, (G/G)o < 4 x 10 -10 yrs u ' 12 ). 

For the Bianchi type / model, by using expression for u(4>) given by 
Eq.flSU) in Eq.flSSD along with values f5 = -2, a > 6/5, 7 = 0, t Q = 14 ± 2 
yrs and m = 1/2, we obtain (G/G) = 1.4287 x 10~ 10 . It also safely lies 
within the allowed range of variation of G for cosmic acceleration. Thus our 
obtained models satisfy the observational limit of G for cosmic acceleration. 



6 Summary and Discussion 

This paper investigates the possibility of obtaining cosmic acceleration by 
using the role of BD parameter in the presence of viscous and barotropic 
fluids. For this purpose, we consider FRW and BI universe models. The 
constructed models entirely depend upon the values of the parameters a, (3, 7 
and m. Firstly, we discuss the FRW model in the presence of viscous fluid. 
We see that the total effective pressure contains a negative factor associated 
with bulk viscosity which leads to negative effective pressure. Consequently, 
the fluid acts as a dark energy candidate and can explain many aspects 
of evolution of the universe. The deceleration parameter constraints the 
parameter a for cosmic acceleration i.e., a > 1. 

For the constant coefficient of bulk viscosity, we obtain (3 = and = 
—2a. The first case leads to GR while for the second choice, in all eras 
of the universe except vacuum dominated era, the BD parameter u(t) is 
a decreasing function of time with small negative values. In the vacuum 
dominated era, we see that for viscosity greater than 0.11, it is possible to 
achieve cosmic acceleration with positive values of u(t). For the variable bulk 
viscosity coefficient with n — 1/2, u(t) corresponds to decreasing function as 
— t 2 with small negative values for different small values of viscosity coefficient 
£0 and — 1 < 7 < 1. This gives rise to accelerated expansion of the universe 
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for a > 1. For the radiative fluid, we have found that the coefficient of 
viscosity appears in exponential function. Here u(t) is a decreasing function 
with negative values both for the accelerated and decelerated phases of the 
universe (—3/2 < uj < 0). 

Secondly, we have taken the BI universe model in the presence of perfect 
fluid with barotropic EoS. Here we have taken two cases when V(<f>) = and 
V((f>) 7^ 0. In the first case, when u(<f)) = Uq, by taking different possible 
choices for the parameters, we see that the BD parameter takes small negative 
as well as positive values for — 1 < 7 < and certain ranges of a. Thus 
the cosmic acceleration can be achieved for positive larger values of to with 
different values of a. Also, for the present universe with lu = —5/3 (taken 
from the negative observed range for cosmic acceleration —2 < u < —3/2), 
we have found <p(t) ~ t~ 3 . In this case, the acceleration rate of the universe is 
higher than Bertolami et al. (<p(t) ~ t~ 2 ) n ) and Benerjee et al. (0(t) ~ t~ 5 / 2 ) 
12 \ When to = to((j)), taking different values of the parameters, we see that for 
all phases of the universe, the BD parameter u{4>) have small negative values 
and lies within the range u < —3/2 which corresponds to cosmic acceleration 
and in agreement with already found results u \ 

For V((f)) 7^ and 00 ((f)) = u>o, we have evaluated the values of scale 
factors A, B, scalar field and V(4>). It is found that in all phases of the 
universe, these values of scale factors A(t) and B(t) lead to q < for all 
positive constant m with m 7^ 1 and a > 1 which corresponds to cosmic 
acceleration. Finally, for V{<p) 7^ and co = w(0), we see that V{<p) is a 
decreasing function starting from zero with the increasing except for the 
choice (5 = — 2 with particular values of 7 and a > 1.3. These values provide 
positive potential energy as they are increasing function of (j). However, for 
the constraint j3 = 2 — a(m + 2)(l + 7), it is possible to have positive potential 
energy for larger values of a in matter dominated era while for smaller values 
of a in radiation dominated era. 

It is worthwhile to mention here that all models discussed here satisfy the 
observational constraints for the variation of Newton's gravitational constant 
available in literature 11 > 12 ) which provides a support to our obtained results. 
Although in each case, we have explained the phenomena of cosmic acceler- 
ation for different ranges of the corresponding parameters. However, these 
ranges of the BD parameter, except for few cases, are incompatible with solar 
system constraints which require u > 40, 000. This is the generic problem 
noted in the context of scalar tensor theories. It would be of great interest 
to see whether this problem can be resolved using other Bianchi models. 
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In order to check the viability of dark energy models based on modified 
theories of gravity, the evolution of cosmological perturbations and the back- 
ground expansion history of the universe may be studied. This can be done 
through Chameleon and Vainshtein mechanisms which suppress the propa- 
gation of fifth force and provide consistency with local gravity experiments 
50,51) _ Q ne ma y ac [ pt these procedures to check the viability of above dis- 
cussed models. 
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